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ADMITTING THE INADMISSIBLE: ADJOINT FORMULATION FOR INCOMPLETE 
COST FUNCTIONALS IN AERODYNAMIC OPTIMIZATION * 

EYAL ARIAN* AND MANUEL D. SALAS* 


Abstract. We derive the adjoint equations for problems in aerodynamic optimization which are improp- 
erly considered as “inadmissible”. For example, a cost functional which depends on the density, rather than 
on the pressure, is considered “inadmissible” for an optimization problem governed by the Euler equations. 
We show that for such problems additional terms should be included in the Lagrangian functional when 
deriving the adjoint equations. These terms are obtained from the restriction of the interior PDE to the 
control surface. Demonstrations of the explicit derivation of the adjoint equations for “inadmissible” cost 
functionals are given for the potential, Euler, and Navier-Stokes equations. 

Key words, adjoint, aerodynamic optimization, Euler, Navier-Stokes, optimal shape. 

Subject classification. Applied and Numerical Mathematics 

1. Introduction. In recent years there has been a growing interest in solving optimization problems 
governed by the Euler and the Navier-Stokes (NS) equations (for example [1]-[10]). The new interest in 
this classical field is due to advances in computer performance and improvements in algorithms for the 
numerical solution of the flow equations. Among the many optimization methods that are being pursued, 
the Lagrange multiplier method or adjoint method is particularly attractive, because of its efficiency for 
problems with many design variables. The adjoint method is based on a variation analysis of the Lagrangian 
and requires that the variation vanishes at the optimum. This necessary condition yields an optimality 
system of coupled PDEs consisting of the state equation, the costate or adjoint equation, the optimality 
condition, and boundary conditions for the state and costate equations. 

Recently, Anderson and Venkatakrishnan [8] reported some difficulties in the derivation of boundary 
conditions for the costate equation for certain cost functionals. The same problem was later reported in [9]- 
[10]. In [8], for example, it is concluded that for aerodynamic optimization problems which are governed by 
the compressible Euler flow equations, only cost functionals that depend solely on the pressure, F = F(p), are 
admissible; and for viscous flow, using the compressible (NS) equations, only cost functionals which involve 
the entire stress tensor (e.g., drag) are admissible. As stated in [8], the difficulty with the “inadmissible” 
cost functionals stems from the need for a suitable balance between the different terms in the variational 
form of the Lagrangian; for some cost functionals such a balance does not exist and the requirement that 
the variation of the Lagrangian vanish does not result in a boundary value problem for the adjoint variables. 
In [9], the authors suggest, for the compressible NS with adiabatic boundary condition on the solid wall, 
to introduce a contribution into the cost functional which depends on the temperature so that appropriate 
cancellation in the variation of the Lagrange functional will occur. In [10], the authors concluded that no 
other choices of cost functionals, other than those suggested in [8], lead to a well-posed problem. 

From the theory of functional analysis costate variables exist for all cost functionals [11]. This, however, 

•This research was supported by the National Aeronautics and Space Administration under NASA Contract No. NAS1- 
97046 while the authors were in residence at the Institute for Computer Applications in Science and Engineering (ICASE), M/S 
403, NASA Langley Research Center, Hampton, VA, 23681-0001. 
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is not the same as saying that all cost functionals lead to a proper boundary value problem for the costate 
equation, which is what we mean by an “admissible” cost functional. 

In this paper, a general method is presented for formulating proper boundary value problems from cost 
functionals considered “inadmissible” in the literature. The method avoids redefining, or introducing new 
terms in, the cost functional. The term “inadmissible” is obviously incorrect. We claim that for so called 
“inadmissible” cost functionals additional auxiliary boundary equations are needed in the Lagrangian. These 
relations are obtained from the restriction of the interior PDE and its derivatives (up to the highest order 
possible) to the boundary. With these additional relations, proper cancellation of terms in the variation of 
the Lagrange functional can be obtained for any well defined cost functional. However, there is value in 
distinguishing those cost functionals that lead to a proper boundary value problem with out the need of 
auxiliary boundary equations from those that do not. To that end, we define complete-cost functionals as 
those that lead to a well posed boundary value problem of the costate equation with out need for augmenting 
the Lagrangian functional with auxiliary boundary equations and we define incomplete- cost functionals as 
those requiring the use of auxiliary boundary equations. 

The paper is organized as follows. In §2 we begin with the potential equation and an admissible cost 
functional. This example is intended to illustrate the problem that arises later with the cost functional 
containing which is also treated in this section. In §3 the adjoint equations are derived for a cost 
functional containing the density, p ) for shape optimization problems governed by the compressible Euler 
equations. In §4 the adjoint equations are derived for a cost functional containing only the pressure, p, 
for shape optimization problems governed by the compressible NS equations. In §5 we discuss our findings 
and make some concluding remarks. Appendix A contains the definition of the Euler Jacobian matrices. 
Appendix B contains identities of polar coordinates which are used extensively in §§3-4. Appendix C contains 
a demonstration of the adjoint derivation on a more complex cost functional, than the one presented in §4, 
for shape optimization problems governed by the compressible NS equations. 


2. The Potential Equation. Let H be a two dimensional domain confined in the area between two 
circles with radii R\ < R 2 > We denote by T the circle with radius R\ and by dfi, - T the circle with radius 
R 2 . Let /*(£) be a given L 2 (F) function defined on the boundary T, and let the function ar(f), defined on 
T, be the design variable. 


2.1. An “Admissible” Cost Functional. Consider first the derivation of the adjoint equation on a 
Dirichlet cost functional using the standard procedure. The minimization problem is defined as follows: 


( 2 . 1 ) 


minFi(0) = 



2 

da 


subject to 


A</> = 0 in 

(2-2) §£ = t on T 

cj) — (f>o on dfl — r 

where ti an d t, respectively, denote the outer normal and tangential directions on the boundary T . The 
definition of the state equation (2.2) assumes the following smoothness properties: 


aeH l ( T) 

4>c Hi (n), 
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where H k denotes the Sobolev space of order k (see [12]). The necessary conditions for a minimum are 
derived with the adjoint method. We introduce the Lagrange multipliers A defined on the domain 0, ( 
defined on part of the boundary T, and d defined on the rest of the boundary 9ft — I\ In terms of the 
Lagrange multipliers, a Lagrangian is defined by: 

(2.4) C(<fi,a,\,C$)= J r (<t>-r) 2 da- + £c(^-^)da + - 4>)da. 

The variation of the Lagrangian is given by: 


(2.5a)£ = j r 24,(4, ~r)d°-j n i&XdSl + f m (* + 1 <( 1 ; " ^) d ” + L-r ***' 


(2.5b) - - iAXdSl + J^(~ + 2(4, - n)*> + J r §|j (C - X)d* + ^ 


The term containing <j> on <9f2 — T was omitted since <j> is fixed on that part of the boundary (see Eq.(2.2)) 
and therefore its variation vanishes. 

The necessary conditions for a minimum are obtained by requiring that the integrands in (2.5b) vanish. 
Matching the terms that multiply <^>, and those that multiply results in the following equations: 


m-. 

AA = 0 in Cl 

m- 

fi = -2 ( 4 >-n onr 

S(T): 

C = A onT 

(dQ - T) : 

A = 0 on dfl — r 

d(r): 

ft= 0 onr. 


Therefore, the costate (adjoint) boundary value problem is defined by 


( 2 . 6 ) 


AA = 0 in Q 
£ = -2 (4>-n onr 
A = 0 on dQ — T, 


with the Frechet derivative of the cost functional with respect to a given by 


(2.7) 


dFi 

da 


dX 

dt 


on r. 


2.2. An “Inadmissible” Cost Functional. Suppose that we want to minimize a cost functional that 
depends on the second normal derivative on the boundary. Then the minimization problem is defined as 
follows: 


( 2 . 8 ) 


minF2 

Q 




2 

do 


subject to Eq.(2.2). However, the definition of the cost functional (2.8) assumes that the second normal 
derivative of <p on the boundary exists and is in L 2 (T). This is not consistent with the smoothness requirement 
of </>, (2.3), unless we additionally assume that a is smoother than required by the state PDE, i.e., 


(2.9) 


aeH 2 ( r) 

4>eHi{ n). 
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If we derive the necessary conditions for a minimum as we did for the Dirichlet cost functional ( 2 . 1 ), then 
the variation of the Lagrangian is given by: 


< 2io) i= L 2 U{U- r ) dc - L ( MA ) dn + L + 1 c ds • f )*'■ 


(As previously discussed, the term that contains (f) on the “far- field” , 5S7 — T, was omitted.) 

In the previous example, we saw how the first and third terms in the variation of the Lagrangian (2.5a) 
combined to give a boundary condition for along T. In this example, the first term in (2.10), the term 
with can not be combined with any of the other terms on the boundary T and, hence, we can not obtain 
a boundary condition on T for the costate equation. A cost functional exhibiting this behavior was termed 
“inadmissible” in the literature [ 8 ]. However, since the state equation ( 2 . 2 ) is linear, <f) depends on a linearly, 
the cost functional F 2 is a quadratic in a. A quadratic cost functional has a unique minimizer, thus the cost 
functional F 2 is in- fact admissible. 

In the next subsection, we show how to overcome this problem using an auxiliary boundary equation. 


2.3. Auxiliary Boundary Equation (ABE). For convenience, let us introduce polar coordinates 
(r, 9). In polar coordinates, the restriction of the interior PDE ( 2 . 2 ) to the boundary T results in the 
following ABE (assuming the consistency requirements of (2.9)) 


( 2 . 11 ) 



on r. 


Here we assume that r — —h is perpendicular to the boundary, pointing into the domain, and R\ is the 
radius of curvature. By Eq.( 2 . 11 ) we get (dt = R\d6) 


( 2 . 12 ) 



on r, 


which can be used to replace in the variation of the Lagrangian. Integration by parts, along the boundary, 
of the term containing tangential derivatives results in a variational form of the Lagrangian which contains 
only the naturally occurring boundary terms (0 and §^). 

Another way to use Eq.(2.11), which leads to the same result, is by adding it to the Lagrangian with a 
new Lagrange multiplier 77 . Then t is augmented with the term /r — -$) da which results 

in the following adjoint equations on the boundary: 

0(0 = -77 + 2(0 -r) = 0 

(2- 13 ) 0(0 : -C?-fr = ° 

-f£(0: ^77-A + C = 0. 

Therefore, the adjoint equations can be written in the following “strong” form: 


(2.14) 


AA = 0 in ft 

~ %£) onr 
A = 0 on 3f2 — r. 


The Frechet derivative of the cost functional with respect to a is then given by 


(2.15) 


dFj _ cK d\ 2 / d 3 <P df*\ 

da dt dt R\ \dr 2 dt dt ) 
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In the same manner other cost functionals can be treated by taking (if necessary and consistent with 
the smoothness assumption of the cost functional definition) higher derivatives of, A0 = 0, restricting the 
resulting equations to the boundary and adding them to the Lagrangian with additional Lagrange multipliers. 

We add a remark on the required smoothness of the Lagrange multipliers. If we insist in solving the 
adjoint equation in its “strong” form (2.14), then the smoothness requirement of the state solution should 
be even stronger than (2.9) (also f*(s) should be smoother than L 2 (F)). This requirement can be relaxed 
if the adjoint equations are solved in a weak formulation, for example, by finite elements [11]. For future 
reference we give here the following definition. 

Definition 2.1. We define a complete-cost functional as one that leads to a boundary value problem 
of the costate equation with out the need for augmenting the Lagrangian functional with auxiliary boundary 
equations . Otherwise , the cost functional will be termed as incomplete . 

In the next sections we treat similarly higher level models of fluid dynamics PDEs. In these cases, the 
functional analysis of the smoothness properties of the state equations are not as obvious as in the above 
example. 

3. The Euler Equations. For simplicity the derivation is done in two dimensions. Let U denote the 
vector of state variables: 


(3.1) 


V = (p, pu, pv, pE ) T , 


where we use the following notation: p denotes the density, u = (u, v) denotes the velocity vector, and E 
denotes the total specific energy. 

The Euler equations are given by (conservative form) 


(3.2) 


div(pu) = 0 
div(pu 0 u + pi) = 0 
div(puH) = 0 


where u ■ n = 0 on the solid wall T, and with additional appropriate boundary conditions on the far-field. 
(The above system, in the interior, can be written in an equivalent form in terms of the Jacobian matrices, 
V(AU) = 0, where A = (A, B ) denote the Jacobian matrices given in appendix A.) 

The following are state relations for the pressure, p ) and the total enthalpy, H: 


(3.3) 


p=(y-l)pE- VpN 2 
H — -3—2 + Ml 

3.1. Natural Boundary Terms. Integrating by parts the Euler equations results in a term containing 
the pressure on the boundary: 


(3.4) 


j V(AU)dQ= f p(0, ni, ri2, 0) T cf<7, 
Jn Jan 


(in polar coordinates n\ = — 1 and — 0), therefore any minimization problem that contains terms other 
than the pressure will result in non-canceling terms in the variational Lagrangian. 

3.2. Example of an Incomplete Cost Functional. The following cost functional is incomplete (the 
design variable is the shape of the solid wall T): 


(3.5) 


F(p)= j^(p-p-) 2 da. 
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The Lagrangian is given by 

(3.6) C = f r (p — p*) 2 dcr + A T V(AU)dQ -f f r C(u • n)da- f far-held terms, 
where A is the vector of costate variables, 

(3.7) A = (Ai , A 2 , A 3 , A 4 ) t 
and we denote by A the adjoint “velocity” vector, i.e., 

(3.8) A = (A 2 , A 3 ). 

The variation of the Lagrangian is given by 

(3.9) C = f r 2 p(p - p*)da + J q U(A t K x + B T A y )d^l + f r p\ 2 dcr + f r ag(U , K)da+ far-held terms, 

which results in non- cancellation of terms in £, since the cost functional is not given in terms of the natural 
boundary term p L 

However, in general we can write 


(3.10) p = p{p,s) 

where s is the entropy and, in particular, on the surface T we can write 


(3.11) 


P- Pi?) 


in the absence of shock waves 2 , since the entropy is constant along the streamline wetting the surface. 
Hence, we know that we can overcome this problem. In the next section we derive the adjoint equations, for 
the incomplete cost functional (3.5), with a general procedure, along the lines of §2.3. 


3.3. Auxiliary Boundary Equations. As in the potential problem we derive auxiliary boundary 
equations (ABEs) by restricting the interior PDEs to the boundary. For simplicity, we examine the result- 
ing equation locally around a point on the boundary and use polar coordinates. Fig.(l) depicts the local 
coordinate system. Throughout the paper we use the unit tangential vector t on the boundary instead of 6. 
Note that on the boundary, T, 


(3.12) 


d_ _ l_d_ 
dt ~ Rd6' 


Also, in polar coordinates the components of the velocity vector u will be denoted by u = (u r ,u*) where 
u T — u ■ t and ut = u • 6\ similarly, the components of the adjoint “velocity” vector A will be denoted by 
A = (A r , A*) where A r = A • t and At = A • 6. 


3.3.1. The Continuity Equation. In polar coordinates, on the boundary, the continuity equation is 
given by (V(pu) = 0 ) 

(3.13) pu r ) + ^(pu t ) = 0 . 

Also higher order derivatives of the continuity equation can be taken and restricted to the boundary and 
are considered auxiliary boundary equations, as long as the solution in the interior is smooth enough so that 
these derivatives exist. 

1 In fact, for shape optimization problems, as we discuss here, the variation of the Lagrangian includes more terms on the 
boundary T that depend on 5 ■ ft. However, these terms contribute only to the gradient term j r ag(U,K)da and therefore do 
not play a role in the derivation of the adjoint boundary value problem. For simplicity we do not discuss these terms in this 
paper. 

2 In the presence of shocks, it is still valid to write p — p(p) in a piecewise sense between shocks using the Rankine-Hugoniot 
conditions to connect the piecewise regions along the streamline wetting the surface. 
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Fig. 1. The local coordinate system around a point on the boundary. 


3.3.2. The Momentum Equations. The momentum equations, V • ( pu 0 u + pi) 
coordinates are given by (see Eq.(B.4)) 


(3.14) 


JU^r) + ^ + \-§e{pUrU t ) ~ ^ ^ - 0 

l(pu r u t )+2 ^ + i^(p^) + f| = 0 . 


0 , in polar 


Using the solid wall boundary condition and the relation (3.13), the restriction of the above to the boundary 
results in 


(3.15) 


pu t 

R 


+ dr 


ou . 4 . §2 — n 
T a t dt^dt u > 


where R denotes the local radius of curvature (see Fig. ( 1 )). 

Higher order equations can be derived by taking derivatives of equations (3.14) and restricting them to 
the boundary. 


3.3.3. The Energy Equation. The energy equation in polar coordinates is given by (V * ( puH ) = 0) 

( 3 .i6) ~{putB) + ~{pnirB) = 0. 

Using the solid wall boundary condition and the relation (3.13), the restriction of the above to the boundary 
results in 


(3.17) 



or 


H — const on T. 


3.3.4. The Derivation of the Adjoint Equations. The definition of the cost functional (3.5) con- 
tains an implicit assumption that the restriction of the density state variable to the solid wall is continuous 
(and also that /?(T) E L 2 (r)). However, in general the density is not continuous in the direction perpendicu- 
lar to streamlines while the pressure, p, and the normal velocity, u r , are always continuous in that direction 
[13]. The new assumption on the smoothness of the density is introduced into the Lagrangian by adding the 
ABEs: 


(3.18) C ABE = £ [m . n) + £(pu ■ i)) + m ( - ^ + g) + m {put ^ + §?) + 

where ( 771 , ■ ■ * , 774 ) are additional Lagrange multipliers. For the cost functional (3.5), the restriction of the 
continuity and first momentum equations are not required and therefore we choose 771 = 772 = 0 . The 
variation of the enlarged Lagrangian, £ + Cabe , yields three adjoint equations on the boundary for A, 773 
and 774 : 


(3.19) 


P( r ) : +2(p-p*) = 0 

Ut(F) : ~ ^(pwt%) + = 0 

p(r) : + -V = 0 
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where we have used the relation 3 

(3.20) (p " + “*** = °- 

The above system can be solved by solving the first two PDEs in (3.19) for 773 and 774 and substituting the 
result in the third equation which is the transpiration boundary condition for A. 

4 . The Navier-Stokes Equations. The compressible NS equations are given by 

div(pu) = 0 

( 4 . 1 ) div(pu 0 u) = div(a) 

div(eu) + div(q) = div(au) 

where the stress tensor, < 7 , is given by 


<7 = — pi + fi 2 div(u)I -f pdef(u), 


and I denotes the unit tensor, def = grad + grad T , and p and p 2 are the first and second viscosities, 
respectively (/i 2 + \p = 0). The vector q denotes the heat conduction vector, q — - kgrad(T ), where T 
denotes the temperature, and k denotes the coefficient of conductivity and will be set equal to a constant. 

I ^jjl 2 np 

The total energy satisfies e — p 

The solid wall boundary conditions are given by 


(4.2) 


u — 0 

aT + b^ 


where a, 6 and c are parameters (in this paper we set a = c = 0 and 6=1, resulting in the adiabatic 
wall boundary condition). 

4.1. Natural Boundary Terms. For simplicity we will denote the system of NS equations by V-F = 0. 
where F consists of the flux vectors. Integration by parts of the compressible NS equations results in the 
following boundary terms: 

(4.3) f (V • P)dfl = f (0, (an)i,{<rn)2,0) T d(T. 

J n Jdn 

Therefore, the natural boundary terms for the compressible NS equations are the total fluid force components, 
(<rn)j. (In other words, the only complete cost functionals are those which measure lift or drag.) 

4.2. Example of an Incomplete Cost Functional. Let us take, for example, the following cost 
functional which is incomplete (here, as in the previous section, the design variable is the shape T) since its 
variation is not given in terms of the force components in (4.3): 

(4.4) F (p) = jjjP-P* ) 2(ia ■ 

The cost functional (4.4) was treated previously in [14] by neglecting a term in the Lagrangian and in [8] by 
modifying the cost functional. A demonstration of the adjoint derivation on a more complex cost functional 
is given in appendix C. 

3 we assume that the term defined on the endpoints of T which results from the integration by parts, J r ^3 |f dir, is equal to 
zero, i.e., r)3P(9T + ) — T)sp(dV~) = 0. 
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The Lagrangian is given by 


(4.5) £ = f r (p — p*) 2 da + f Q A T ■ (V * F)d£l + / r C^da + / r d • ud<7+ far-field terms. 


The variation of the Lagrangian is given by 
£ = interior terms + f r 2 p(p — p*)dcr+ 


(4.6) 


fan ~Tk^ + %k\ 4 + + f*/xA t + pA r da + / r d 5 (C/, A )do. 


The term ^ does not contribute to the adjoint equations (assuming adiabatic boundary conditions) and 
therefore will be omitted (by the same reasoning the term containing the variation of u on T and the variation 
of the far-field terms vanish). Also note that the variation in T can be transformed to variations in p and 
p using the equation of state, T — The cos ^ functional (4.4) is considered incomplete since the term 

determines the adjoint boundary condition A r = 0. As a result we do not obtain a boundary value 
problem representation of the adjoint equations. 


4.3. Auxiliary Boundary Equations. 


4.3.1. The Continuity Equation. Since the boundary condition on T implies that §i{put) = 0, we 
get from the continuity equation that (see Eq.(3.13)) 


(4.7) 



4.3.2. The Momentum Equations. Let us write the momentum equations in the form 


div(pu 0 u) = div(cr) = div(-pl) + div(r). 

The term div(pu 0 u + pi) is given in polar coordinates in Eq.(3.14). Taking the limit to the boundary we 
get that 


lim div(pu 0 u) = 0. 


The term r in polar coordinates is given by 


(4.8) 


= p, 2 div(u)I 4* pdef(u ) = 


H 2 div{u) + + 


+ ^div(u) + 2 M (^ + ^) 

Evaluating div(r) on the boundary and using the identities 

limf__r div(u) = 0 


(4.9) 


1 dp du* 


fb div (v) 

results in the following auxiliary momentum boundary equations (see Eq.(B.4)) 


(4.10) 


d 2 Uj 
dr 2 




J -n^± — (I 
R^ 1 dr ~ U * 
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4.3.3. The Energy Equation. The energy equation is given by 


div(eu) + div(q) = div(au). 


Using the relation 


eu 




ur 

P^r + 




2 ’ 7(7-1)' 

and taking its divergence in polar coordinates we get 
(4.11) 

The term div(q) = —kAT satisfies (assuming an adiabatic boundary condition) 

<d 2 T d 2 Ty 


lim div(eu) = 0. 


(4.12) 


Umdiv(f) - -*(^ 3 - + w )- 


The term au and its tangential derivatives on the boundary are zero, therefore only its radial derivative 
is considered: 

(du t \ 2 


(4.13) 


liin div(cru) = p 


fouty 

\dr) ’ 


Eqs.(4.11)-(4.13) imply the following auxiliary energy boundary equation 


(4.14) 


d 2 T Q2 T fOutY 
k \dr 2 + dt 2 ) Hdr ) 


0. 


4.3.4. The Derivation of the Adjoint Equations. As in the Euler case, we add the ABEs to the 
Lagrangian: 


f r (du r \ (dp 1 dpdut 
Cabe = J r h ( aT ) + ~ 


dp 


1 dp du t 


d 2 u r d 2 u t 
dr 2 




drdt ) 


(4.15) 


(dp d 2 u t 1 dut\ ( d*T\ (cfu t \ z \1 , 

+7?3 (at " r^W) + n - Ha? + a?) “ K )r fT * 


<9 2 T d 2 T> 


/9u t \ 2 N 


However, for the cost functional (4.4) we need only the first term in Eq.(4.15) and therefore we choose 
7/2 = m = r?4 = 0: 

C + Cabe = interior terms 4- j r 2p(p — p*)da+ 

(4.16) [-fk^ + (»7i + f/x-V) + + pA r j da + J r agi(U, A, rn)da. 

The variation of the enlarged Lagrangian, £ + Cabe, yields the following equations on the boundary (using 

T=^- l (-jiP+ 1 p p)) 


^(r): 


Vl + 


(4.17) 


m- 2(p-p*) + A r -^ T ^=0 


fKr) = 

m- 


p\t = o 
fc 1 9A* 

7 — 1 ~p* dr 


0. 


The resulting adjoint equations on the solid wall are given by (See Eq.(4.2)) 
(4.18) 


A r = — 2(p - P*) 
At = 0 


d\ 4 
cfn 


= 0. 
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Note that the first equation in (4.17) determines the auxiliary Lagrange multiplier 771: 

4 8 

(4.19) r)i = — -/zA r - -p(p - p*). 

That Lagrange multiplier is affecting the gradient, g\ = g\ (U, A, 771), since the variation of the continuity 
auxiliary boundary equation (4.7) with respect to a change in the shape, T, results in variational boundary 
terms that multiply only the term a. 

5. Discussion and Concluding Remarks. We present a method for the derivation of the costate 
equations for problems in which the cost functional does not lead to a proper boundary value problem for 
the costate equation, when derived in the standard way. We define such cost functionals as incomplete; 
it is required to “complete” the Lagrangian with auxiliary boundary equations in these cases in order to 
derive a boundary value problem of the costate equation. We demonstrate the method on three problems 
involving the incomplete cost functionals using the potential, the compressible Euler, and the compressible 
NS equations. Our aim is to give the costate equations a representation of a boundary value problem and 
not to treat rigorously the issue of existence of solutions to the resulting system of costate equations. We 
note that for all cost functionals it is possible to derive the adjoint equations in the discrete level in the 
standard manner; the problem of incompleteness exists only in the PDE level. The relation between the 
costate equations that we derive and the discretely derived costate equations in the limit of mesh-size going 
to zero is beyond the scope of this paper. 


Acknowledgments. The authors thank Dr. Bambang I. Soemarwoto, of the Institute of Technology 
Bandung in Indonesia, for fruitful discussions on the derivation of the adjoint equations for the Navier-Stokes 
equations. 
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(A. 1 ) 


Appendix A. Definition of the Euler Jacobian Matrices. 

/ 0 1 

(3 - 7 )u 

—uv v 


-u 2 +^(u 2 + v 2 ) 


0 0 \ 

— (7 — l)f 7 — 1 
u 0 


-u( 7 £-(7-1)(u 2 + ^ 2 )) 7 E-*^(3u 2 + v 2 ) -(7-1 )uv 7 « ) 


(A.2) B = 


0 

— uv 

L zl 
2 


-i; 2 + ^{u 2 + v 2 ) 


0 

v 

-(7- 1 )u 


1 

u 

(3 - i)v 


0 \ 

0 

7 - 1 


-v(yE - (7 - l)(u 2 + v 2 )j — (7 - l)uu jE - ^(u 2 + 3v 2 ) 'yv j 

Appendix B. Identities of Polar Coordinates. 

(B- 1 ) W(/) =(§£>; %) 


(B.2) 


du r 1 du t u r 


(B3) 


grad(u) — 


du r J_ ® u t 

dr r dd 

dxtt 1 _j_ u. 




dr 


dd 


If A = Aij is a tensor then 


(B.4) 


divA = 


1 All _t_ 1 dA i2 

dr ' r ' r d& 


dA 


dAu 

dr 


1 A21 1 I dA 2 2 1 Al 

' r ‘ r d6 ' r 
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Appendix C. Demonstration of the Adjoint Derivation on a Complex Cost Functional 
Governed by the Compressible Navier-Stokes Equations. 

Let us take, for example, the following cost functional which is incomplete 

(ai) F = / r<? -p-)^ + / r (|-r)V 

The Lagrangian is given by 

(C.2 )£ = f r [p - p da + / r — f* j da + f Q K T ■ (V ■ F)dQ + / r + / r far-field terms. 

The variation of the Lagrangian with respect to a change in the shape, d, is given by 
L — interior terms + f r 2 p^p - p*Jda + f r 2§£ — f* jda+ 


(C- 3 ) 


f F rrn ^4 9T 8u r 4 dHt . _ 1 , 

L l ~ sT + ^ 1 + aTS^ + + H‘ i7 ' 


For the cost functional (C.l) we need only the first two terms in Eq.( 4 . 15 ) and therefore we choose 773 = 
rj 4 = 0: 


(C. 4 ) l ABm [u (tpO + *(§!- 


_ 1 <9p 1 dp du t 1 dp dut d 2 u r 

p 2 dt dr + p dt dr + p dt dr ^ dr 2 




drdt) 


The term doesn’t cancel with any other term in the Lagrangian and, therefore, we add another ABE 
by taking the radial derivative of the continuity equation, 


(C. 5 ) 

and restricting it to the boundary: 
(C.6) 


d 2 d 2 

a^ (pUr) + a^ (pUt) = 0 ’ 


9 d 2 u T dpdut , d 2 u t n 

l,m -d „ <*.) = p -^- + -— +p —= 0 . 


We can further add to the variation of the Lagrangian (Eq.(C. 3 ) added with Eq.(C. 4 )) the variation of the 
residuals of (C.6) with a Lagrange multiplier 775: 


(C. 7 ) 


£new = £ + CaBE + 


IA 


d 2 U r 


+ P 


d 2 u r dpdu t dpdut „d 2 u t d 2 u t 


dr 2 ' r dr 2 dt dr dt dr ^ drdt ' r drdt 

The resulting adjoint equations on the boundary are given by 4 (using T = ~-j- ^ — ^p + pPj) 


+ P-, 


;) 


da. 


^(r): 

^(T): 

f(r): 

P(F): 


7 Jl + %p\ r = 0 
— 2/U772 + T] 5 p = 0 

\ k 1 dA 4 r\ 

Ar ~~p dr “ U 


pAf + 772M2^|f + - £ §j? i = 0 

p( r ) : 2 (p-p*) + -7?2M24rff^ - P2 J 7 (r72^^) +P5f^ = 0 . 

The above system can be solved by first solving for 771 , 772 and 775 from the first three equations and then 
substituting the result in the last three equations, which are the desired adjoint boundary conditions. 


4 we assume that the terms defined on the endpoints of T which result from the integration by parts are equal to zero. 
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